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Abstract 

Including corrections of order 0{rriK* /mg), we present an analysis of pho- 
tonic penguin contributions to the decay B —* K*j in the perturbative QCD 
framework. Employing several models of the meson wave functions, we demon- 
strate that the corrections of 0(rriK* /tub) are enhanced and will provide sub- 
stantial contributions to the decay because of the B meson wave function being 
sharply peaked (bound state effect). The numerical predictions for the correc- 
tions are about 30% ~ 60% which depend on the non-perturbative inputs such 
as the meson wave functions and the 6-quark mass. 
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1. Introduction 

The rare decay B — > X*7 has attracted great attentions especially after the CLEO 
Collaboration first identified this decay and gave its branching ratio PJ. The decay 
B — > i^*7 is dominated by the flavor- changing quark-level process b — > S7 which 
can occur not only through penguin diagram at one-loop level in the standard model 
(SM) but also through virtual particle in the supersymmetry and other extensions of 
the standard model p], |3|]. Thus accurate experimental measurements and theoretical 
calculations of this decay can provide a precision test of the standard model as well 
as a test of new physics at present experimentally accessible energy scale. It has been 
pointed out [|J that perturbative QCD (PQCD) may be applicable to the exclusive 
nonleptonic decays of B meson since there is a hard-gluon exchange between the heavy 
and light quarks in these decays. Recently calculations also show that PQCD may give 
a good description of the two body hadronic decays of B meson || . 

In the standard model (SM), the mainly contribution to the decay B — > if* 7 comes 
from the photonic penguin diagrams which are shown in Fig. 1. Compared to Fig. la, 
the contribution from Fig. lb is of order m^./mg, and thereby it is not included in 
Ref. ||. In this paper we shall re-analyse decay B —* K*^ in the SM by including 
0(m,K* /tub) corrections in the amplitude. Including bound state effect and employing 
several models of distribution amplitudes of B and K* mesons, we find that these 
corrections are enhanced by the bound state effect and become more important. This 
paper is organized as follows: In section 2, we calculate the photonic penguin diagram 
contributions to the order m^./ms in the amplitude after describing the effective 
Hamiltonian. In section 3, we present numerical results by employing several models 
of meson distribution amplitudes. As usual, the last section is reserved for summary. 

2. Contribution coming from photonic penguin diagram 

The effective Hamiltonian (the square blob part in Fig. 1) which describes the 
photonic penguin diagram, can be expressed as |7], |§L |S| 

Heff = -^V tb V*CM0 7 (ii), (1) 
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where 

°M = -^- 2 m bS^ u F,J-{l + 75)6. (2) 

In the above expressions, CV(/z) is the Wilson coefficient which contains the effects of 
QCD corrections, 



cm = v -^o \c 7 (m w ) - i| (V 0/% - 1) - ^ (v 28m - 1) 



(3) 



where r\ = a s (fi)/a s (mw), Po = H — (2/3)?t./ and Cj{mw) = —0.19 is given in the 
W-mass scale. 

The wave function of the B meson can be written in the form 0] 

1>b = ^^=0b(>) 7s(i^ - m B ), (4) 

where I c is the identity in the color space. For the K* meson, the wave function can 
be expressed as 

i>K* = ^ -J^>k< (x)f (jb K * + m K * ) , (5) 

where £* is the polarization vector of the K* meson. 0g and 4>k* are the distribution 
amplitudes of the B and K* meson respectively. 

We express the contribution to the amplitude in the gauge invariant form 2 

Mi=Ux -— 
2ps ■ q 

The contributions from Figs, la and lb can be written as 

1 1 



Pb ■ qt* -e-PB- tq ■ C + ie^cplfetfe* ?* • (6) 



ti = Gm b [dx][dy](j) B (x)(f) K *(y) 



1 ~/2 7 71 ° F ^ ^ h h)lc 

AGm b [ [dx] — 4> B (x) [ [dy] ^ 1 V ^ m f — 2mbl7lB ^ ( y ) 
Jo xi Jo y x [ml - (1 - yi)m B \ 



2 It is worthwhile to note that the expression for the amplitude presented here, Eq. (g), is gauge 
invariant, while the one given in El is not because of the second term in the bracket of Eq. (Q) being 
missed ]To} ] . 
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t An f u l 1 a f\ f u l m 6-2(l -yi)m B }m K * 

+AGm b / [dx\—0 B {x) / [dy\ ^ <Pk*KV) 

Jo xi Jo yi mi - 1 - vi mi 



xi jo j/i [mg - {i - yijm^J 
4Gm b I 1B I™ + 4Gm b I 1B I? K L ° 

LO , + iVLO 



and 



t^o = Gmb l[dx][dy]<j> B (x)<l> K .(y) ' 



/ 6 2 - m\ k 2 g 



isijfa - m B ) J*U>k* + m K .) , y y ta , F 1 
^2 \71 (fb + m b j(T p -r^2 ( 



-AGm b [ [dx] p-(j) B {x) / [dy] — </>*-*(?/) 

-4Gm 6 J 2 B/ 2 ^ , 



In the above expressions, [dx] = <ixidx2<5(l — xi — X2), [dy] = dyidy 2 S(l — yi — y 2 ) , 9 
and e are the momentum and polarization of the photon respectively, and 

G = ^£-V tb V*C F C 7 (fi)ea s (fi). (9) 

Xi and yi in Eqs. (0) and (|J) are the momentum fractions carried by the light quarks 
in the B and K* mesons respectively. The distribution amplitude of B meson, <p B (x), 
should be sharply peaked at some small value of X\ since m b is much larger than the 
light quark mass [Q, [|. Thus we keep only the leading contributions of xi in the quark 
and gluon propagators, which are X\ terms in Fig. la and x\ terms in Fig. lb. The 
fermion propagators in Figs, la and lb contribute different factors to t% and t 2 : The 
one in Fig. la involving only K* meson variable in the form of \ j\y\m\ — (m B — ml)} 
is attributed to the integrals l[%* and The one in Fig. lb involving only B 

meson variable in the form of l/(xim B ) is attributed to the integral l 2B . The gluon 
propagators in Figs, la and lb involving both B and K* meson variables in the form 
of l/(xiyim B ) can be factored to the integrals Ii B and 1^. In this way, t\ is factorized 
to two independent integrals U B and Uk*- 

In Eqs. (0) and t\° provides leading contribution while t^ LO and t^ LO are cor- 
rections of 0{m K * /m B ). It is interested to notice that the suppression factor m K */m B 
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in t^ 10 can be compensated by the bound state effect as it is going to be demon- 
strated in the following. Compared to Iib, the fermion propagator in Fig. lb provides 
an additional factor \jx\ to Iib- Because the distribution amplitude of B meson, (f>B, 
is sharply peaked at x± ~ 0.05 ~ 0.1 f|, I<ib is much larger than I\b- For example, 
employing a simple model for (fi B , <Pb ~ ~ e_e) with 



the ratio is (see Table 1) 

This factor will cancels approximately the suppression factor mjf./mg being about 
1/17 in I2K* 1 which make the contribution coming from Fig. lb become important. 
There is no similar enhancement factor in t± , so it is order mp/mB and may be 
neglected as compared to t\° . 

It has been pointed out || that the contribution coming from Fig. la, contains a 
large imaginary part because of the pole of the heavy quark propagator. This imaginary 
part does not correspond to the long-distance physics. It should be noticed that the 
ratio of the imaginary part to the real part depends on the 6-quark mass m b (namely 
€b) and B(K*) distribution amplitudes (see Table 2), which are about 3.5 ~ 0.8. Thus 
the contribution form Fig. lb should be taken into account although it provides only 
a real contribution. 

The decay width and branching ratio can be obtained readily, 

r = ^ 

167rmB 

Br(B - iT 7 ) = (13) 

J- total 

3. Numerical calculation and model analysis 



m B -m b 

£b = , (10) 



mb «1Q~18. (11) 
es rn B - m b 



E ( M l + M 2 



polarization 



(12) 
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For the numerical results, we take the following parameters as inputs: 



tb 



A Q cd = 200 MeV, 
m w = 81 GeV, 
V tb = 0.999, 
/ B = 132 MeV[ll], 
1.46 x 10~ 12 second. 



fi = 1 GeV, 
m t = 2mw, 
V u = -0.045, 
f K . = 151 MeV0, 



(14) 



The numerical results should depend on the expressions of distribution amplitudes 
4>b(%) and 4>K*{y) which are determined by the non-perturbative physics. For the <pB 
we adopt the following models: i) According to Brodsky-Huang-Lepage prescription 
the B meson wave function can be given in the form flT4" |, 



iPb{x, k±) = Aexp 



_ b 2 ( m 2 b + k\ t m 2 g + k\ 



+ 



x 2 Xx 

in which the parameters A and b are determined by two constraints: 



(15) 



(16) 



and 



B 



\- dx A^ \^b(x, k ± )\ 2 SB 
o Ibn 6 



(17) 



Pb is the probability of finding the \qq) Fock state in the B meson. The second 
constraint Pb ~ 1 is reasonable since with the increase of the constitute quark mass the 
valence Fock state occupies the most fraction in the hadron, and in the nonrelativistic 
limit the probability of finding the valence Fock state is going to approach unity. Then 
we can obtain the distribution amplitude of B meson 



A 



;X\X2exp 



_ h 2(< + ™£ 



lQlT 2 b 2 \ 22 X\ 

ii) Szczepaniak, Henley and Brodsky suggested another model for (Pb(x) 0, 

.4 



(18) 



(e|/xi + l/x 2 - 1) 
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(19) 



where A and e# are given by Eqs. (|T^) and (|T0|) respectively, iii) The simplest model 
for 4> B is the 5-function approximation which has been adopted in Refs. || || 

<j> B {x) = ^=8{x x -e B ), (20) 

where eg is related to the longitudinal momentum fraction of the light quark (see Eq. 
(0))- 

We adopt the following two models for (fix*'- i) it has been pointed out |14[ |l|, [Uj 
that K* meson wave function is close to its asymptotic behavior, so we adopt the 



expression in Ref. |14| 



V 2/2 Vi , 



(21) 



where A = 41.4 GeV -1 , b = 0.74 GeV _1 , m s = 0.55 GeV and m g = 0.35 GeV. The 
quark masses appearing in the meson wave functions (distribution amplitudes) should 
be the constituent quark masses since the wave function is determined mainly by the 
soft-physics, while the quark masses appearing in the hard amplitudes should be the 
current quark masses which can be ignored reasonably for the light quarks, ii) The 
asymptotic expression for 4>k*, 

<pK*(y) = V3f K * yi y 2 . (22) 

The numerical results are given in tables 2 and 3. <Pb HL i <Pb HB an d 4>b nave 
different behavior in the x-space (see Fig. 2). is not so sharply peaked as <p B HB 

and 4> S B , and the position of the maximum of 4> B HL is farther away the end-point x± — 
than that of the other two models i.e. does not emphasize the small-xi region 

so strongly as (j) B HB and <p B . Thus the value of I\b [J-2b) calculated with (pB HL is the 
smallest one among the three models (see table 1). Because of ti and t 2 depending 
on <^ B only through the integrals Iib and I 2 b respectively (see Eqs. ([7p and (^)), the 
decay amplitude and branching ratio calculated with 4>b HL will be also the smallest 
one (see table 3). 

It can be found that t^ LO is about 1/10 of t LO because of the suppression factor 
rriK*/mB, while t^ 10 is the same order as the real part of t\ since the bound state 
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effect compensates approximately the suppression factor m K */m B (see table 2). The 
corrections to the decay branching ratio are about 30% ~ 60% which varies with the 
distribution amplitudes of B and K* mesons and the 6-quark mass. The corrections 
calculated with <f) B HB is more important than that with the other two models, and 
the corrections calculated with and (j)™* are very similar since they have similar 
behavior. It can been found also that the corrections become more important with 
m h increasing. We would like to point out again that it is because the distribution 
amplitude of B meson should be sharply peaked at some small value of X\ = e B (the 
bound state effect) that the corrections of 0(m#* /ms) coming from Fig. lb become 
more important. 

As comparing with the experimental data, we find that the results calculated with 
4> B HB and (f> B , and with rrif, being about 4.9 are comparable to the experiment data 
Br(B -> K*j) = 4.5 ± 1.5 ± 0.9 x 10~ 5 0. 

4. Summary 

The decay B — > K*^ is a very attractive process since it provides an experimentally 
accessible way for a subtle test of the standard model and a test of new physics. Both 
more accurate theoretical calculations and more accurate experimental measurements 
about this decay mode are worthwhile and necessary. By including the corrections of 
order ms-./mg in the photonic penguin diagrams, we have analysed the decay B — > 
K*^ in the framework of perturbative QCD. Employing several models of the meson 
wave functions, we find that the 0(m#* /mg) corrections coming from Fig. lb provides 
substantial corrections to the branching ratio since the bound state effect provides 
an enhancement factor mj/(mg — mj) which cancels approximately the suppression 
factor m^*jm B . The corrections are about about 30% ~ 60% which depend on the 
non-perturbative inputs such as the meson wave functions and the 6-quark mass. 
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Foundation of China. 
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4> B B HL 


<t> S B HB 


<t>B 


m 6 (GeV) 


4.8 


4.9 


5.0 


4.8 


4.9 


5.0 


4.8 


4.9 


5.0 


hB 


0.29 


0.30 


0.31 


0.38 


0.46 


0.59 


0.42 


0.53 


0.72 


I2B 


3.50 


3.63 


3.77 


7.16 


10.7 


18.4 


4.20 


6.83 


12.8 



Table 1. (^-dependence of U. 







<t>% 


m 6 (GeV) 


4.8 


4.9 


5.0 


4.8 


4.9 


5.0 


<P B b HL 


+ LO 
l \ 


-1.40-2.291 


-1.87-2.291 


-2.44-2.121 


-0.68-2.321 


-0.97-2.441 


-1.35-2.661 


+ NLO 
L l 


0.01-0.291 


-0.06-0.311 


-0.14-0.301 


0.12-0.281 


0.08-0.321 


0.03-0.381 


-iNLO 
L 2 


-0.79 


-0.84 


-0.89 


-0.86 


-0.91 


-0.97 


<P S b HB 


+LO 
L l 


-1.80-2.941 


-2.83-3.451 


-4.63-4.031 


-0.88-2.881 


-1.47-3.691 


-2.58-5.051 


+ NLO 
L l 


0.01-0.371 


-0.08-0.471 


-0.27-0.581 


0.15-0.361 


0.13-0.501 


0.06-0.721 


+ NLO 
L 2 


-1.63 


2.47 


-4.33 


-1.76 


-2.69 


-8.13 


<& 


+LO 
L l 


-1.96-3.241 


-3.24-3.961 


-5.63-4.871 


-0.96-3.151 


-1.68-4.241 


-3.13-6.151 


^NLO 


0.01-0.411 


-0.10-0.531 


-0.33-0.701 


0.16-0.401 


0.15-0.561 


0.80-0.871 


+ NLO 
l 2 


-0.95 


-1.57 


-3.02 


-1.03 


-1.72 


-3.29 



Table 2. Decay amplitudes in units of 10 8 GeV. 
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<t>V* 


m 6 (GeV) 


4.8 


4.9 


5.0 


4.8 


4.9 


5.0 


<P B b HL 


Br LO 


0.60 


0.73 


0.87 


0.46 


0.58 


0.74 




0.95 


1.20 


1.50 


0.70 


0.91 


1.21 


g r Full _ g r LO 


37% 


39% 


42% 


34% 


36% 


39% 


<t> B HB 


Br LO 


0.99 


1.66 


3.16 


0.76 


1.32 


2.69 


Br Full 


1.89 


3.70 


8.90 


1.40 


2.83 


7.16 


g r Full _ g r LO 


48% 


53% 


65% 


46% 


53% 


62% 


4>B 


Br LO 


1.19 


2.19 


4.65 


0.90 


1.74 


3.96 


Br Full 


1.80 


3.71 


9.35 


1.33 


2.82 


7.47 


g r Full _Q r LO 
B r Full 


34% 


41% 


50% 


32% 


38% 


50% 



Table 3. Branching ratio Br(B — ► K*j) in units of xlO 5 . 
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Figure Caption 

Fig. 1. The photonic penguin diagram. The square blob represents the effective 
vertex. 

Fig. 2. The distribution amplitudes of B meson employed in our calculation: 0^ HL (the 
solid curve) with m h = 4.9 GeV and m q = 0.35 GeV; (fig HB (the dashed curve) 
with e B = 0.072; <j) s B ~ 5(x l - 0.072) is not plotted in this figure. 
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Fig. 1 
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